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Interaction between a galactic disk and a live dark halo 
with an anisotropic velocity distribution 
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Abstract. We have extended previous analytical studies of the interaction of dark halos with galactic disks by 
introducing for the halo particles anisotropic distribution functions in phase space. For this purpose we have 
employed the shearing sheet model of a patch of a galactic disk embedded in a homogeneous halo. We find that 
velocity anisotropy increases considerably the maximum growth factor of perturbations in the disk. 
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1. Introduction 

The dynamical evolution of non-axisymmetric structures 
such as spirals or bars is driven by the redistribution of 
angular momentum within the disk galaxy that harbours 
them. This was argued and demonstrated by Lynden-Bell 
& Kalnajs (1972), who showed that stars at inner Lindblad 
resonance and at higher order resonances within corota- 
tion emit angular momentum, which is absorbed by stars 
at corotation, at outer Lindblad resonance and at higher 
order resonances outside corotation. Non-resonant stars 
can also contribute to this exchange, provided the spi- 
ral or bar perturbation has a non-negligible growth rate. 
Since the non-axisymmetric structures are inside the coro- 
tation circle negative angular momentum 'perturbations' 
(Kalnajs 1971, Lynden-Bell & Kalnajs 1972), this ex- 
change induces their growth. In this early picture the role 
of the spheroid is not included. However, recent, state-of- 
the-art A'^-body simulations have clearly shown that it can 
be predominant in galaxies with strong bars. Indeed, if the 
galactic disk is embedded in a spheroid or dark halo, a sim- 
ilar angular momentum exchange between the bar and the 
surrounding bulge or dark halo is expected. Athanassoula 
(2002, 2003) studied in particular the interaction of a bar 
in the disk and a live spheroidal component (bulge and/or 
halo) . Provided the latter has a distribution function that 
is a function of the energy only, it is possible to show ana- 
lytically that it can at all its resonances only absorb angu- 
lar momentum. Thus in the presence of spheroid or a halo 
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a bar can grow stronger than in its absence, since in the 
former case more angular momentum is taken from the 
bar region. Tremaine & Weinberg (1984) and Weinberg 
(1985) calculated that the angular momentum exchange 
also causes a considerable slow-down of the bar, which 
is similarly enhanced in the presence of a live spheroid 
with an isotropic distribution function. This was, at least 
qualitatively, confirmed by a number of A'^-body simula- 
tions (e.g. Little & Carlberg 1991, Hernquist & Weinberg 
1992, Athanassoula 1996, Debattista & Sellwood 2000, 
Athanassoula 2003, O'NeiU & Dubinski 2003, Valenzuela 
& Klypin 2003). Athanassoula (2002, 2003) used A-body 
simulations to show that a considerable amount of halo 
material is at resonance with the bar. She further traced 
the angular momentum exchange in the simulations and 
found very good qualitative agreement with analytical 
predictions. Angular momentum is emitted by particles 
within corotation: mainly by particles at inner Lindblad 
resonance, but also by particles at higher order resonances 
and also by non-resonant particles. It is then absorbed by 
disk particles at corotation and outer Lindblad resonance, 
halo resonant particles at all resonances and non-resonant 
halo particles. Indeed, bars in simulations with a respon- 
sive halo, that can take part in the angular momentum 
exchange, grow considerably stronger than bars in simu- 
lations with a halo which has the same mass distribution, 
but is rigid (Athanassoula 2002). This justifies, in retro- 
spect, the doubts voiced much earlier by Toomre (1977) 
as to the entirely passive role of the halo in bar formation 
and development. 



2 



B. Fuchs and E. Athanassoula: Anisotropic live dark halos 



The dynamics of bars and spiral density waves in galac- 
tic disks are closely related to each other (Toomre 1981, 
Bertin at al. 1989a, b). Accordingly enhanced growth due 
to the presence of a responsive halo is also known for spi- 
ral density waves. Mark (1976b) has studied the angular 
momentum transfer from mode like spiral density waves 
to a bulge - halo system surrounding the galactic disk. 
The density waves were modelled as rigidly rotating expo- 
nentially growing spiral patterns whose dynamics was de- 
scribed in the l/FiCi?J approximation (Mark 1974, 1976a). 
The models presented in Mark (1976b, hereafter referred 
to as M76) show considerably enhanced growth rates of 
the density waves due to angular momentum loss to the 
bulge - halo system. Recently Fuchs (2004, hereafter re- 
ferred to as FG4) followed a different approach to the spiral 
amplification question in the presence of a live dark halo. 
Using the shearing sheet model he showed how a live dark 
halo responds to swing amplified spiral density waves and 
found also a very enhanced growth of such density waves. 
These density waves interact with the halo via halo par- 
ticles either on orbits in resonance with the waves or on 
non-resonant orbits. 

The shearing sheet model (Goldreich & Lynden-Bell 
1965, Julian & Toomre 1966) has been developed as a tool 
to study the dynamics of a patch of a galactic disk. This 
is assumed to be infinitesimally thin and its radial size is 
assumed to be much smaller than that of the disk. Polar 
coordinates can be therefore rectified to pseudo-Cartesian 
coordinates and the velocity field of the differential rota- 
tion of the disk can be approximated by a linc^ar shc^ar 
flow. These simplifications allow an analytical treatment 
of the problem, which helps to clarify the underlying phys- 
ical processes operating in the disk. 

Up to now all studies of the effect of a live dark halo on 
the dynamics of non-axisymmetric structures in a galac- 
tic disk have assumed an isotropic velocity distribution 
of the halo particles. On the other hand, it is well known 
from cosmological simulations of dark halo formation that 
the velocity distribution of the halo particles must be 
anisotropic. If the anisotropy is measured by Binney's pa- 
rameter (3 = 1 — {(Tg + CT^)/2crr (cf. Binney & Tremainc 
1987) this ranges from /3 « near the centers of the halos 
to /3 « 0.5 at the virial radii (Colin ct al. 2000, Fukushige 
& Makino 2001). Hansen & Moore (2005) have suggested 
that the radial variation of /3 is roughly linear to the log- 
arithmic density gradient of the halo, (5 cx —dlnp/dlnr. 
Very recently Abadi et al. (2005) have reported a simula- 
tion of a mildly anisotropic dark halo with an anisotropy 
parameter of /3 « 0.3. In the same experiment they have 
also simulated the formation of the baryonic galaxy in- 
side the dark halo. Interestingly, the velocity anisotropy 
of the spheroidal component of the model galaxy is signif- 
icantly larger, (3 « 0.4 to 0.8, which is consistent with the 
anisotropy parameter of /? 0.7 of Galactic halo stars in 
the solar neighbourhood (Arifyanto et al. 2005). Following 
a different concept, which describes dark halos as par- 
tially relaxed dynamical systems, Trenti & Bertin (2005) 
find also radially anisotropic velocity distributions with 



an anisotropy parameter in the region around corotation 
of /3 « 0.6. Helmi et al. (2002) discuss the implications 
of such an anisotropy for experimental searches of dark 
matter particles. 

Thus the studies of the interaction of dark halos with 
non-axisymmetric structures in galactic disks have to be 
extended to halo models with anisotropic velocity distri- 
butions. As a first step we consider in this short note the 
amplification of spiral density waves in disks embedded in 
anisotropic dark halos. 

2. Introducing anisotropic distribution functions 
for the halo particles 

2.1. Response of a live dark halo to a density wave in 
the shearing sheet 

We study now the response of a live dark halo with an 
anisotropic velocity distribution to a swing amplified den- 
sity wave in the disk. For this we use the method and nota- 
tion of F04, assuming that the shearing sheet is immersed 
in a homogeneous live dark halo. Since we modify only the 
unperturbed distribution function of the halo particles in 
phase space, we can immediately follow the formalism of 
F04 to derive the Fourier transformed Boltzmann equa- 
tion which describes the evolution of the perturbation of 
the distribution function /hk (eqn. 20 of F04). We adopt 
now the anisotropic background distribution function 



/hO = 



Ph 



2n cr^a^ 



exp 



2a 



(1) 



u and V denote the planar velocity components of the halo 
particles with v pointing in the direction of the differential 
shear in the sliect. w is the velocity component perpen- 
dicular to the plane of the sheet. As in F04 we choose 
orthogonal Cartesian coordinates ^, r], and C with ^ paral- 
lel to the wave vector k of /hk- According to the symme- 
try of the distribution function (1), one of the coordinates 
perpendicular to k can be aligned without loss of general- 
ity with the V velocity component. Thus the distribution 
function (1) can be written as 



/hO = 



xexp 



Ph 



27r' cr^fTv 



■ exp 



v"^ (I'll cos a — tij, sin a) 



{v\ I sin a + v± cos aY 
2^E 



(2) 



fc2 



1.2 



where cos a = k\\/ + A:^, and k\\ = |k||| = 

is the component of the wave vector k parallel to the mid- 
plane of the sheet. v± denotes the second velocity compo- 
nent perpendicular to k. The Boltzmann equation can be 
then immediately integrated with respect to v and, after 
some algebra, the integration with respect to leads to a 
Boltzmann equation which has the same form as eqn. (20) 
of F04, but where the one-dimensional velocity dispersion 
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ah of the isotropic distribution function is replaced by the 
velocity dispersion of the velocity component in the ^- 
direction parallel to k 



Ceff = 



COS + (t£ sin 



(3) 



We can proceed then exactly as in F04 and obtain the 
same final formal result $hk|| (z = 0) = T(A:x, fcy, a;)<i>dk|| , 
where the function T is modified by replacing cth by a^n 
in eqns. (15) and (17) of F04. 

As alternatives to the distribution function of the halo 
particles (1) we have considered also the anisotropic dis- 
tribution functions 



hO 



and 



Pb 

3 ~ ^'^P 



277 a^al. 



Pb 



(4) 



(5) 



If these are integrated with respect to the two velocity 
components perpendicular to k, we find the effective ve- 
locity dispersions 



CTeff = \/ 0-2 cos /32 -|- cr^ sin /32 



and 



(6) 



(7) 



respectively. The angles /3 and 7 are given by cos/3 = 



fcj^ . An effective ve- 



fcx/ ^k'^ + A:^ and cos 7 = fcy/ ^Jkf 

locity dispersion can also be calculated analytically 
for a triaxial distribution function. However, the result is 
so cumbersome that we do not consider this case here. 



fT||(a|| -FcTw) fc2 + fc;^ 

+4:TT^Gph exp 



(8) 



t-t' 



2Ak'y 



I fcx fcx 



:e^^[-\^al-al){kl + k'l)(t-t' +^-^) ] 



xerfc( ^^kl + k'/ 



2Akl 



The modified equation (35) of F04 can be integrated nu- 
merically still with very modest numerical effort. The 
characteristic response of the shearing sheet embedded in 
a live dark halo to an initial impulse is not changed qual- 
itatively by the anisotropy of the velocity distribution of 
the halo particles. The shearing sheet develops swing am- 
plified density waves with their amplitudes enhanced by 
the presence of a live dark halo as compared to a static 
halo (cf. Fig. 2 of F04). In the case of the distribution 
function (4) we obtain the Fourier transform 



/ 



a^kl + a^k'lia^^k^ + ^'y + y ^u^x + <^w^'y) 
,~ / / k-v /t„ 



2Ak' 



-4:TT^Gp\y exp 



/ , k^ — k'^ 



t-t 



2Akl 



(9) 



2.2. Dynamics of the shearing sheet embedded in an 
anisotropic live dark halo 

The halo response to the perturbation in the disk has to 
be inserted into the Boltzmann equation describing the 
evolution of the distribution function of the disk particles 
in phase space. Wc follow again the formalism of F04 and 
find the fundamental Volterra type integral equation for 
the Fourier transforms of the perturbations of the gravi- 
tational potential of the shearing sheet in the same form 
as eqn. (33) of F04, but with a modified Fourier trans- 



form J-(Te ^""-y )t-t'. The latter can be calculated an- 
alytically as in F04 using first formulae 6.317 and 3.952 
of Gradshtcyn & Ryzhik (2000) and then formulae 3.223 
and 3.466 for the integrals with respect to fcz, leading to 



j T(fcx,fcy,w)e* 



:exp[-^(a^-a^)fc^{i-t' + ^) ] 



xerfc(^V^^ + fcr 



2Ak'^ 



and similarly for the distribution function (5) 



T T(fcx,fcy,a;)e' 



STT^Gpb 



^alkl + (Tlk'l{a^ ^kl + k'l + ^alk^ + a^k'j) 
kv /ijv 



xSlt-t' + 



2Ak^ 
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2Ak' 



(10) 
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xexp[-l(„J-ai)t'j((-(' + %5^)'] 

xerfc(^V«^h'' + ^|)- 

The characteristic response of the shearing sheet embed- 
ded in dark halos with these anisotropic distribution func- 
tions of the halo particles to an initial impulse is still the 
same as in the previous case, only the maximum growth 
factors of the density waves change. All three Fourier 
transforms (8,9, and 10) are reduced to the form given in 
F04, if the distribution functions become isotropic. 

2.3. Numerical applications 

2.3.1. Comparison with previous work 

No results with anisotropic halos, either analytic or nu- 
merical, have been published so far. Indeed, as mentioned 
already in the introduction, the lack of such information 
has been one of the main motivations for our work. A cou- 
ple of previous works, however, have studied the influence 
of a live dark halo with an isotropic distribution function 
on density waves and we will present here some possible, 
albeit rough, comparisons. 

The interaction of a WKBJ type density wave with 
a separate non-rotating component was studied for the 
first time by Marochnik & Suchkov (1969). The astro- 
nomical community was at that time not aware of the ex- 
istence of dark halos, but the considerations of Marochnik 
& Suchkov (1969) carry over to the disk - dark halo inter- 
action. Unfortunately, only very few concrete details were 
given. Much more detailed is the investigation by Mark 
(M76). In Fig. 3 of M76 it is shown how much stronger 
the amplification of the density waves is in the presence of 
a live halo than in an isolated disk. This depends on the as- 
sumed corotation radii of the spiral patterns. Since swing 
amplified density waves are basically corotating with the 
center of the shearing sheet, we adopt for each corota- 
tion radius of the spiral modes the parameters given in 
Table 1 of M76 for that galactocentric radius and calculate 
for these parameters the corresponding maximum growth 
factors of swing amplified density waves in the shearing 
sheet. The regions in the galaxy models, where the ro- 
tation curves are essentially flat, are most suitable for a 
comparison. Using the formalism discussed in the previ- 
ous sections, we find for the parameters of model Ai at R 
= 6 kpc an enhancement of the maximum growth factor 
of the amplitudes of swing amplified density waves due to 
the halo by a factor of 2, while Fig. 3 of M76 indicates 
that the growth rates of the density waves are enhanced 
by a factor of 2.3. Similarly, we find for the parameters of 
model Bi in the range R = 14 to 18 kpc maximum growth 
factors enhanced by factors of 3.2 to 1.4. Fig. 3 of M76 in- 
dicates for that range an enhancement of the growth rates 
by a factor of 1.7. Thus the numerical values are in rough 
agreement. 




0.3 0.4 0.5 0.6 0.7 



ij 

Fig. 1. Maximum growth factor of the amplitudes of 
density waves in the shearing sheet as function of the 
anisotropy of the velocity distribution of the halo parti- 
cles. The solid line with open circles illustrates the case 
if the components of the velocity dispersion are varied as 
c^u = cTy These are drawn as abscissa. The third com- 
ponent is chosen as a^f = -^/l — 2a^, so that the to- 
tal velocity dispersion is constant. The cases = CTw, 
(Tu = -^/l — 2(7^ are shown as the dashed line with open 
squares and = cry,, = y'l — 2a^ as the dotted line 
with filled asterisks, respectively. The isotropic case a = 
0.58 reproduces the result of F04. 



2.3.2. Numerical applications to anisotropic halos 

We have determined maximum growth factors of the am- 
plitudes of the density waves for various values of the 
parameters characterizing the disk - halo configuration. 
Representative results for the three cases described above 
are presented in Fig. 1. In all three cases we have taken 
<^tot = CTu + CTv + CTw = 1> = 0.5 and k'^ = 

0.5 fccrit- In F04 it was shown how equations (8) to (10) 
can be cast into dimensionlcss form. The velocity disper- 
sion of the halo particles is then measured in terms of the 
velocity dispersion of the disk stars. We have assumed for 
the examples in Fig. 1 a ratio of dtot '■ ^disk = 5:1. The 
velocity dispersion of the disk stars is determined by the 
Toomre stability parameter for which we have adopted a 
value of Q = 1.41. For the density of the dark halo and the 
circular velocity of the center of the shearing sheet around 
the galactic center we have adopted values typical for the 
solar neighbourhood in the Milky Way, namely p^G/k'^ = 
0.01 and vq^o '■ Cdisk = 220 : 44. We find in all three cases 
the same characteristic behaviour, that the amplification 
of density waves in the shearing sheet is smallest, if the dis- 
tribution function of the halo particles is almost, but not 
exactly, isotropic. Even in the isotropic case, though, the 
amplification is still considerably larger than in a shearing 
sheet embedded in a static halo. The dependence of the 
maximum growth factor on the velocity dispersion com- 
ponents is rather subtle. If we concentrate on the case 
(Tu — <Jv, it becomes clear from equation (8) that rises, 
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if cr|| is decreasing with respect to CTw, which leads to a 
stronger input to the swing amplification mechanism. This 
can be intuitively understood as the halo becomes more 
susceptible to the sheared trailing spiral arms in the disk, 
if the planar velocity dispersions become small. The de- 
pendence of T on the CTw velocity dispersion component is 
twofold. Technically, this is related to the integrations of 
the halo response to the disk perturbations with respect 
to the fcz wave number (cf. eqns. 26 and 28 of F04). The 
result is that T rises also at lower dw velocity dispersions. 
This means that the dark halo becomes then less stiff in 
the vertical direction so that it can support the density 
waves in the disk more effectively, although the effect is 
less pronounced. The relative minimum of amplification 
of density waves is at slightly larger (7^ = velocity dis- 
persions than the isotropic case a^ — a^ — a^ — 0.58 fXtot- 
The other diagrams in Fig. 1 can be explained by similar 
arguments. However, since in these cases only one planar 
velocity dispersion is varied, the maximum growth fac- 
tor is less than in the previous case. We have also con- 
sidered other numerical values for the parameters Q and 
A of the disk as well as for the halo parameters pb and 
(7tot • Even though the overall amplification of the density 
waves can change drastically, the effect of the anisotropy 
as compared to the isotopic case remains always roughly 
the same. 

Extrapolating these results, we predict that spirals 
and bars growing in disk galaxies with non-isotropic ha- 
los should be on average stronger than those in isotropic 
halos. Unfortunately it is not possible to test this for 
real galaxies. It is, however, possible to do so for A^- 
body simulations and such a study is presently underway 
(Athanassoula, in prep.). 
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